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Recent progresses in the understanding of the scaling behavior of self-avoiding flex- 
ible polymerized membranes (tethered manifolds) are reviewed (joint works with 
B. Duplantier and E. Guitter, and with K. Wiese). They rely on a new general 
renormalization group approach for a class of models with non-local singular in- 
teractions. This approach allows to prove the existence of a e-expansion for the 
scaling exponents, and validates the one loop results obtained by direct renormal- 
ization methods. Applications of the method to polymerized membranes at the 
tricritical 0-point are presented. 



1 What are polymerized membranes ? 

Polymerized membranes, also called tethered membranes, are an interesting 
model for two-dimensional flexible polymerized films □. In particular, they 
generalize the well-known problem of the statistical behavior of long chains, 
i.e. of polymers. However, even simple properties, such as the statistical 
properties at thermal equilibrium of a single membrane, are quite difficult to 
study theoretically, and require the development of new techniques. L shall 
review some recent progresses made with B. Duplantier and E. Guitter □, and 
some applications of these techniques with K. Wiese □. 

Let me consider an ideal flexible two dimensional network fluctuating in 
three dimensional space, characterized by: (i) a fixed and regular internal 
structure with a fixed connectivity (for instance a triangular lattice), and (ii) 
negligible bending rigidity. The only forces between elements of the network 
are the elastic forces (described for instance by a harmonic potential between 
nearest neighbors). If self-avoiding interactions are not taken into account, one 
deals with a "phantom" membrane, which is known to be, for vanishing or small 
enough bending rigidity, in a crumpled phase. The large scale fluctuations are 
described by a Gaussian model. Such a model is physically different from the 
random surfaces used in models of fluid two-dimensional films (for instance 
for vesicles or micro-emulsions). For these models, the bending rigidity is 
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important and the fluid character of the film is reproduced by a fluctuating 
connectivity of the lattice. 

Taking into account the self-avoiding interactions (which are short ranged 
in the three dimensional embedding space, but involve elements of the lattice 
which can be arbitrarily far apart in term of the internal lattice distance) one 
expects that the long range properties of the membrane will be changed. For 
instance, the average size (R) of the membrane in physical space is expected 
to scale with the internal size L of the lattice as 

(R) oc V < v < 1 (1) 

if the membrane is still crumpled, but swollen by self-avoidance. The mem- 
brane may even become flat (this is observed in numerical simulations of self- 
avoiding tethered membranes in three dimensions); the exponent v is then 
equal to 1. 



2 An Edwards Model for Z?-dimensional polymerized membranes: 

In the continuous model introduced inB@, the field r(x) describes the embed- 
ding of a D-dimensional space, whose points are labeled by the coordinate x, 
into the external <i-dimensional space, whose points are labeled by the vector 
r. The Hamitonian H (the free energy for a configuration r) is the sum of a 
Gaussian elastic energy and of a 2-body repulsive interaction, proportional to 
the coupling constant b: 

H[r\ = J dPx ±(V x r) 2 + b J d D X J d D y 5 d (r(x) - r(y)) . (2) 

The internal dimension D may be taken as a continuous parameter, interpo- 
lating between polymers (D — 1) and membranes (D — 2). The issue is to 
compute with this model critical exponents describing the scaling behavior of 
large membranes, for instance the exponent v (related to the fractal dimension 
df of the membrane by v = D/df), and the configuration exponent 7, related 
to the scaling of the partition function Z of a finite membrane with internal 
extent L by 

Z oc L 7 " 1 constant 1 '" . (3) 

The mean field exponents are obtained by setting 6 = and one obtains 
v = (2 - D)/2 and 70 = 1 - d{2 - D)/2. 

Dimensional analysis shows that mean field theory is invalid if e, the engi- 
neering dimension of &, is positive 

[b] = e = 2D-d(2-D)/2 > 0. (4) 
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Figure 1: d-D plane 



The general picture of the expected behavior as a function of the internal 
dimension D and of the external dimension d is presented in Fig. [I] 

A natural idea is to compute the corrections to mean field by a e-expansion 
a la Wilson-Fisher. Using the so-called direct renormalization, inspired from 
polymer theory B, the first explicit calculations have been performed to first 
order in e by Aronowitz & LubenskyQ and by Kardar & Nelson □ . The ba- 
sic idea of this method is to perform explicit perturbative calculations for a 
finite membrane and for e > 0. Perturbation theory is then UV and IR fi- 
nite, but has UV divergences when e — ► 0. These poles in 1/e can be removed 
by reexpressing the observables in terms of adequate dimensionless renormal- 
ized quantities, such as the 2nd virial coefficient. The internal size L of the 
membrane plays the role of a renormalization mass -1 scale, and the renormal- 
ization group equations can thus be obtained, by writing the L dependence 
of the renormalized theory. The consistency of this procedure at first .order 
has been checked by explicit calculations by Duplantier, Hwa & Kardarcl, but 
its consistency at all orders cannot be proved by the de Gennes trick valid for 
polymers (the scaling limit of the self-avoiding walk can be mapped into a local 
field theory with 0(n) symmetry in external <i-dimensional space, in the limit 
n — > 0, and standard renormalization group theory is then appliquable) . 
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3 Renormalization Theory and Renormalization Group for non- 
local Theories: 



Instead of using direct renormalization, we succeeded in proving the renormal- 
izability of_thc model (Q), by considering it as a non-local field theory in D 
dimensionsa. Perturbation theory is obtained by expanding the observables as 
power series in b. The bi-local "interaction vertex" (in field theoretic language 
I denote it a bilocal operator) is written in Fourier transform as 

• • = 8 d (r(x) - f(y)) = [ d d k e ^(r{x)-r(y)) ( 5 ) 



x y 

It can be viewed in a Coulomb gas representation as the integral over the 
"charge" k of a neutral "dipole" with charge +k at x and charge — k at y. 
The term of order b K in the perturbative expansion of the partition function 
(as well as of other observables) involves K dipoles (xi,yi), . . . {xk,Vk)- The 
integration over the charges k\ , . . . kx gives an integral over the positions of 
the dipoles of the determinant of the "dipole energy" quadratic form Q 



K -d/2 
Y[ dPxidPyi det 
i=i 



Q[xi,yi] ■ (6) 



K 

Q is a K x K matrix such that kiQijkj is the Coulomb energy (in D- 
dimensions) of the K dipoles. Each Qij is a linear combination 

Qij = G (xi,Xj) + G {ui,yj) - G (xi,Vj) - G (x j ,y i ) . (7) 
of the Coulomb potentials Go between endpoints of the dipoles i and j 

G (x,x') = (r(x)r(x')) = \x - xf~ D . (8) 

The Coulomb potential is properly defined for < D < 2 by analytic con- 
tinuation. It is then negative, it vanishes for x — x' (for D > 2 it diverges), 
but as for D > 2 is decreases at large distances. The integration over the 
2K points in non-integer D-dimensional space can also be defined properly by 
analytic continuation in D and the use of distance geometry. This amounts to 
replace the integration over the 2K x D coordinates by an integration over the 
K x (2K — 1) scalar distances between these points. 

One can show that short distance UV singularities occur in the integrals, 
when the determinant det[Q] vanishes. This occurs if and only if some of the 
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Figure 2: A non-neutral regular configuration and two neutral UV singular configurations 



end-points of (not necessarily the same) dipoles coincide, forming "atoms" , and 
so that the dipoles form "molecules, and moreover if one can assign non-zero 
charges fc, to the dipoles while the atoms stay neutral. This condition is more 
easily depicted graphically on Figure 2. 

The associated singularities of these integrals are related to the behavior 
at short distance of the expectation value (with respects to the free Gaussian 
model) of products of bilocal operators as given by Equ. |^. One can show 
that this short distance behavior is encoded in a multilocal product expansion 
(MOPE), which generalizes Wilson's operator product expansion. Let me give 
two examples: 

When the two points x and y of the bi-local interaction operator tend 
towards a single point, this operator can be expanded in terms of local oper- 
ators involving derivatives of the field f. The first terms of the expansion are 
explicitly (omitting unimportant coefficients) 



= \*-yr 2D i 

~x"y +\x- y\ e - D - 2 (x a - y a )(x? - /) : V a r V f : 

+••• (9) 

1 is the identity operator (its expectation value is 1), the : : in the operator 
: V Q r V/3r : denotes the normal ordering subtraction prescription required to 
deal properly with the UV singularities contained in V Q r V^f. 

The second example is less simple, and shows that when the end-points of 
two bilocal operators tend pairwise towards two different points, this generates 
again bilocal operators 



x, 



'.1 — fU ™ |2-D i |„. „, \2-D-\- d / 2 



JZ; 1 = [ki-^| 2 - D + |2/i-y 2 | - — - 

X 2 *' *' y 2 + ... (10) 
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Figure 3: UV divergent configurations and the associated relevant operators 

This structure is generic, and products of local and bilocal operators gen- 
erate multilocal operators of the general form 

r p 

<S>{ Xl ,---,x P } = d% H [(V ?0 p 5 d {r -r{ Xi )) A t { Xi )] . (11) 
J »=i 

where the Ai{x{) are local operators, which can be decomposed into products 
of multiple x-derivatives of r. The m, are integers. For P = 1 and m = 
one recovers local operators A(x) (m > gives 0). For P = 2, mi = mi = 
and A\ — A2 — \ one recovers the bilocal interaction operator, etc. . . These 
operators have a very special form: they can be viewed as a local convolution 
in the external rf-dimensional r space of a non-local product (in the internal 
D-dimensional space) of the P local operators Ai. 

The MOPE implies that the formalism of renormalization theory and of 
renormalization group equations, which has been developed for local quantum 
field theories, can be adapted for this model. One is in fact interested in the IR 
scaling behavior of the lattice model, when some length scale L goes to 00. This 
lattice model is described by the Hamiltonian (0), with a short distance lattice 
cut-uff a. To study this IR limit it is equivalent to look at the UV continuum 
limit of the model when the physical length scale L is kept fixed, while the 
UV cut-off a goes to 0. In this limit one can construct, via renormalization, 
a finite renormalized theory with a — 0, which obeys renormalization group 
equations. From these equations, one recovers the large distance behavior of 
the lattice model I started from. The procedure works well in perturbation 
theory when one is close to the upper critical dimension, i.e. for e small, and 
it leads to the e-expansion. 

In our case, the MOPE can be used to determine, by power counting, 
which multilocal operators are relevant and give UV singularities (poles in 
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Figure 4: Examples of nested singular configurations 



l/e). Then one can also show that these poles can be subtracted by adding 
to the Hamiltonian (|^) counterterms proportional to the marginally relevant 
multilocal operators, leading to the UV finite renormalized theory. For the 
model of self-avoiding surfaces, this analysis shows that the UV divergences 
are associated only with local and bilocal operators, as depicted on Fig. |[ 
and that only three operators are relevant: the identity operator 1, the elastic 
energy operator (Vr) 2 and the bilocal operator S d (r(x) — r(y)). 1 is strongly 
relevant, and gives power-like UV divergences proportional to a~ D (a being 
a short-distance cut-off). The two other operators are superficially relevant, 
they give logarithmic UV divergences or equivalently poles in l/e at e = 0. 

The fact that the so-called superficial divergences, associated to a global 
contraction of points towards a singular configuration, can be subtracted by 
counterterms is a consequence of the MOPE. A complete proof of the renor- 
malizability of the theory is possible, but much more delicate. It requires a 
control of the subdivergences coming from successive contractions associated 
to nested singular configurations, such as those depicted of Fig. || . 

4 Scaling for infinite self-avoiding membranes 

A first application of this formalism is the derivation of scaling laws. Since 
the model is renormalizable (at least perturbatively), it can be made UV finite 
(for e ~ 0) by introducing two counterterms in the Hamiltonian. The new 
renormalized Hamiltonian is of the form 



bn is the dimensionless renormalized coupling constant (the perturbative ex- 
pansion in bft is UV finite order by order). Z is a "wave-function" renor- 
malization factor and Zi, a coupling constant renormalization factor, both are 



H[r\ 




d D x {Vr) 2 + b B f i e Z b 



IS 



d D xd D y S(f{x)-f{y)) . (12) 
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Figure 5: The /3-function and the RG flow for e > 0. 



perturbative series in bu, with poles up to degree l/e K 1 at order K. [i is the 
renormalization momentum scale. As for ordinary local theories, such as the 



Landau-Ginzburg- Wilson <£> 4 Hamiltonian, one can change 6r and r in Equ. 12 
into bare quantities in order to rewrite the renormalized Hamitonian as a bare 
Hamiltonian given by Equ. ^. The renormalization group /3-function and the 
anomalous dimension 7 of the field r are now defined in the standard way 



d 

(3(b R ) = V-Q^ b R 



(13) 



The /3-function is found to be of the form 

0(bn) = —ebn + cbR + 0{b 2 R ) ; c positive constant , (14) 

and therefore there is, at least for small e > 0, an IR attractive fixed point b* = 
0(e), which governs the scaling behavior of self-avoiding polymerized surfaces 
at large distance. The existence of this fixed point ensures the universality 
of this non-trivial scaling for e > Oi, and that no new interactions, possibly 
non-local in external space, are generated by the RG transformations. 

The explicit calculation leads to the same results for the scaling exponents 
v and 7 than the direct renormalization method at first order in e. With this 
method higher order calculations are in principle feasible, but technically quite 
difficult. In particular, already at second order the RG functions cannot be ex- 
pressed analytically, and numerical integration methods have to be developed. 
Work is in progress to compute the scaling exponents at order e 2 . 
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5 Finite Size Scaling and Direct Renormalization 

The model given by Equ. |^ describes an infinite membrane with flat internal 
metric, corresponding to an infinite and regular flexible lattice. Finite mem- 
branes are described by a similar model, but the D-dimensional membrane 
M is now embodied with a fixed non-trivial Ricmannian metric g a p(x), with 
curvature R (examples are closed membranes with the topology of the sphere 
Sd or the torus T D ), and may have a boundary dM (open membrane with 
the topology of the disk for instance). A similar analysis can be performed for 
such models, and the MOPE structure of short distance singularities is still 
valid, but new local operators A{x), which depend on the internal metric on 
M and on the boundary dM, such as the curvature R, appear in the MOPE 
and in Equ. (pi]). The renormalized Hamiltonian now contains at least five 
operators and five independent renormalization factors Z 

H[r\ = [ Zl + [ Z{Vrf + 11 Zb5 d {r-r) 
Jm Jm J Jm 

+ [ ZR + [ Zl . (15) 
Jm JdM 

The curvature operator J M R is superficially relevant only for D — 2 and the 
boundary operator J gM 1 only for D = 1 . When these additional terms are not 
relevant, the first three renormalization factors Z are the same for finite curved 
membranes than for the infinite flat membrane. This property is analogous to 
the renormalization property of local field theories in finite geometries, which 
justifies the finite scaling laws for critical systems in finite geometries, and it 
has two very important consequences: (i) The scaling hypothesis at the basis 
of the direct renormalization approach, which relies explicitly on calculations 
with finite membranes, is shown to be valid to all orders in perturbation theory; 
(ii) for "abstract" membranes with dimension D < 2, with the only exception 
of open polymers (D = 1 open surface), the following hyperscaling relation □ 
relating the configuration and the v exponents holds: 

7 = 1 - vd . (16) 

6 Self-avoiding Membrane at the tricritical O-point: 

Finally, let me briefly discuss recent results obtained with K. Wiese on the scal- 
ing behavior of polymerized membranes at the O-point a. This point separates 
the swollen phase, where the self-avoidance repulsive forces that I considered 
previously dominate, from the dense collapsed phase, where short ranged at- 
tractive forces dominates. At the O-point the effective two body repulsive 
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Figure 6: Relevant interactions for the ©-point in the d-D plane. 



coupling b vanishes, and two different interactions may become relevant. The 
first one is the 3-body contact repulsion, which is usually considered for poly- 
mers 



d D x d D y d D z 5 d (f{x) - r{y)) 8 d (f{x) - f(z)) . (17) 



Z 

The second one is a modified 2-body interaction, repulsive at short range but 
attractive at larger range (A? is the d-dimensional Laplacian) 

• H • = - j j d D xd D y A f S d (r(x)-r(y)) . (18) 

Calculations at first order are not feasible analytically, and already require 
numerical evaluations of complicated integrals. The results of such one loop 
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calculations are schematically depicted on Fig. where the domains where 
the 3-body and modified 2-body terms are respectively relevant are shown. 
This indicates that the last modified 2-body term is the relevant one for 2- 
dimensional membranes in any external dimension d. There is also a quite 
interesting and non-trivial crossover between the two terms around D = 4/3 
d = 6, which must be studied by a double e-expansion. 



7 Conclusion: 

The theoretical study of the scaling behavior of polymerized flexible mem- 
branes leads to the development of new multilocal continuum field theories, 
and to new applications of renormalization group methods. I hope that these 
methods will lead to a quantitative progress in the understanding of the be- 
havior of real 2-dimensional polymerized membranes. This requires results 
beyond first order in the e-expansion (recall that D = 2 correspond to e = 4), 
and a better understanding of the relation between this RG approach and 
more heuristic or approximate methods, such as variational methods or ap- 
proximate recursion relations. The sophisticated renormalization theory for 
multilocal models should hopefully also find applications in other problems of 
statistical physics, or in other areas of theoretical physics. 
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